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ABSTRACT 


This study is concerned with the development of an analytical 
code for Dynamic Stress Analysis and Fatigue Life Prediction of 
Turbomachinery Blades. Emerging blade technologies are finding it 
increasingly essential to correlate blade vibrations to blade fatigue 
in order to assess the residual life of existing blading and for the 
development of newer blade designs. 


On the basis 

of previously done available work 

an attempt 

has 

been made 

to 

refine the approach 

for stress 

analysis 

by 

incorporating 

the 

crucia 1 ly important 

aspect of 

non 1 ineari ty 

in 


damping. The numerical procedure developed for estimation of stresses 
under nonlinear damping is presented. 

A life prediction algorithm, based on the strain-life and 
fracture mechanics concepts to account for the blade fatigue 
phenomenon involving crack initiation and crack growth, has been 
developed as an improvement upon the conventional stress based 
approach. 

The stress analysis code coupled with the life algorithm forms 


an overall package 
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CHAPTER 1 


INTRODUCTION 


Vibration induced -fatigue is a problem of major concern in 
turbomachinery blading. With ever increasing demand for high speed 
turbines, blade vibration and fatigue has become an important subject 
of study. Blade behaviour is strongly influenced by the mean and 
alternating stresses, experienced during operation, by the elastic 
and fatigue properties of the blade material and the chemistry of the 
operating environment. The critical aspect of the problem constitutes 
occurence of large dynamic stresses, due to resonant or near resonant 
vibration at an integral order of engine speed, when a component of 
the vibrating blade matches in both time and space a corresponding 
pattern in the gas/steam stream caused by the nozzle flow field. 

Blade fatigue usually initiates in a region of high stress at 
some structural and/or metallurgical discontinuity. The earlier 
practice had been to avoid fatigue failure by tuning the blade to 
operate away from the natural frequencies. However, in a typical 
turbomachine this may not always be possible, for it may have three 
or four rotors, each having several stages and the turbine may 
altogether have many blades with different characteristics. Dynamic 
stress analysis of the blade, therefore assumes importance. 

Fig. 1.1 summarises the steps involved in the analysis of a 


turbine blade. 
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Fig. 1.1 STEPS IN THE DESIGN OF A TURBINE BLADE 











3 


The major hurdles faced in approaches to blade fatigue analyses 
have been C 1 1 

► a precise evaluation of nozzle excitation forces. 

► obtaining valid damping coefficients for the blade. 

► determination of dynamic stresses under nonlinear damping. 

► analysis of blade behaviour during transient operations like 
start-up and shut-down. 

► definition of a failure prediction model along with evaluation 
of blade fatigue strength properties under broad conditions of 
operating environment. 

1.1 PRESENT WORK 

This study concerns itself with the following two aspects of 
blade design 

► development of a numerical procedure to determine dynamic 
stresses under nonlinear damping. ( Existing data, in 
literature, for typical blades has been employed for 
illustration ). 

► development of a failure prediction model based on strain-life 
and fracture mechanics concepts. 

A software package has been developed which requires the inputs 

• blade geometric and material data 

• excitation data 

• damping data 

• operational data 

• fatigue parameters 


to yield outputs of 

► natural frequencies and mode shapes 

► dynamic stress and displacement fields under forced vibration 
due to nozzle excitation 

► fatigue life estimate 

Starbase graphics support has been used to depict the results 
graphically, wherever necessary. 

Relevant literature has been surveyed in Chapter 2. 

Chapter 3 outlines the various steps in the dynamic stress 
analysis of blades and describes the procedure to incorporate 
nonlinear damping in the analysis. Results are presented for a 
typical turbine blade using the data available in literature. 

The fatigue life estimation model, developed on strain-life and 
fracture mechanics concepts, is given in Chapter 4. Implementation of 
the model is illustrated for assumed cases of blade defects. 

Chapter 5 concludes the present work, outlining the scope for 


future work 



CHAPTER 2 


LITERATURE REVIEW 


Blade vibration has been an area of intensive research. A single 
•free standing blade can be considered as a pre-twisted cantilever 
beam with an asymmetric aerofoil cross-section mounted at a stagger 
angle on a rotating disc. The blade executes vibrations in coupled 
bending-torsion modes, as the centre of flexure may not coincide 
with the centroid for the aerofoil cross-section and coupled 
bending-bending modes due to pre-twist. The problem is complicated by 
other effects such as shear deflections, rotary inertia, warping of 
the cross-section, root fixing and Coriolis acceleration . Various 
researchers have derived solutions to the problem by considering 
individual aspects such as taper, pre-twist, asymmetry of 
cross-section and centrifugal forces, and have made simplified 
assumptions towards secondary effects. They have used two distinct 
approaches towards blade vibration problem : 

• Continuum model approach 

• Discrete model approach 

In the first approach, the energy terms are set up using the 
beam or plate theory ( depending on the aspect ratio ) and the 
solutions to the equations of motion are obtained directly by energy 
methods like Ray leigh-Ritz , Galerkin, Lagrange, Reissner and 
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Nemat-Nassar based on the Hamilton's principle. Other methods of 
solving differential equations of motion of cantilever blades are 
Integral equation approach, Perturbation and Collocation procedures. 
The literature in this area has been reviewed by Rao C23 and 
reference can be made to this paper for a detailed description of the 
contribution made by various researchers. 

Discrete model approaches like Finite Difference Method and 
Finite Element Method are adopted by some researchers. Finite 
Difference Method has been adopted by Carnegie et al C3, 41 for 
rotating asymmetric blades. Sisto and Chang C51 have adopted the 
finite element method for rotating beam type blades. Rawtani and 
Dokainish C6D used FEM to determine natural modes of pre-twisted 


cantilever plates. 

Bossack 

and Zienkiewicz 

C73 extended 

this 

to 

rotating blades. 

Leissa, 

Mac bain 

and Kielb C83 have 

made 

a 

comprehensive study 

of the 

numerous 

previous 

investigations 

on 

free 


vibrations of pre-twisted cantilever plates. 

Blade damping is a complex problem and the main damping 
mechanisms known to influence damping are 

(1) internal friction damping at the root 

(2) material hysteresis damping 

(3) steam/gas dynamic damping. 

The superposition of these mechanisms determines the stress levels 
expected in resonant vibration and establishes the susceptibility of 
the blade for fatigue failure. 

The possible contributions from these mechanisms can be infinitely 
variable and no predictive theory presently exists to relate these 
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damping mechanisms to predict the specific blade damping parameters. 
Lazan C93 has given a comprehensive coverage of material damping. 
Hammons C103 has considered damping arising out of oscillator motion 
of blade in steam environment. Rieger and Beck CL13 have developed a 
rig to simulate centrifugal forces by thermal means in order to 
determine effects of rotation on damping in blades. Rao, Gupta and 
Vyas C 123 designed and built a spin rig with nozzle passing 
excitation simulated by electromagnets. They obtained damping ratio 
as a function of rotational speed and tip displacements for the first 
four modes. Rao and Vyas C133 subsequently developed modal damping 
envelopes for different critical speeds as a function of tip 
displacement. 

The area of forced vibrations has been researched into much less 
when compared to free vibration analysis. Rieger and Nowak C143 have 
used ANSYS program to determine dynamic stresses of a blade group due 
to flow path excitation. Hoyniak and Fleeter C153 made an energy 
balance between unsteady and aerodynamic work and energy dissipated 
due to aerodynamic damping to predict blade resonant vibrations, 
contribution by some researchers in this area. 

Blade fatigue is a problem currently undergoing investigations. 
Rieger C163 has outlined various aspects involved in undertaking life 
analysis of turbomachine blading. Rao and Vyas Cl 3 have used Bagci's 
fatigue surface and Miner's Cumulative Damage rule for a given 
loading history to estimate fatigue life of blades. Rieger C173 has 
proposed the usage of a fracture mechanics approach to determine low 
cycle fatigue life, employing the necessary stress results from the 
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FEM analysis of Rieger and Nowak C 1 4, 183. 



CHAPTER 3 


STRESS ANALYSIS 


The governing blade vibration equations are presented in this 
chapter. Reissner's dynamic -functional in conjunction with the Ritz 
process and Modal Analysis has been employed to -formulate the problem 
o-f -forced vibration o-f the practical case of a tapered, twisted, 
asymmetric, aerofoil cross-section turbine blade mounted at a stagger 
angle on a disc rotating with constant angular velocity <Fig 3.1). 

These equations, as derived by Vyas C193, ignore higher order 
effects such as shear deflection, rotary inertia, warping and 
Coriolis forces. The blade profile data, excitation and damping data 
are taken from the case studies of Vyas C193 for the purpose of 
i 1 lustration . 

The numerical procedure, developed during the course of the 
present study, to estimate blade response under nonlinear damping, is 
then presented. Strength of Materials formulae are employed to obtain 
the stress fields. Results are obtained for stresses as functions of 
various vibrational and operational parameters. 

3.1 EQUATIONS OF NOTION 

Referring to Fig. 3.1, the kinetic energy of the blade has been 


obtained as [193 







11 


■* p \ 


ACCk + r Q) + ( y - r 

y 


L 

©> 2 > dz + i j r e z dz 
x 2 J eg 


l 

pv? C-Jf <x' + 


+ r' ©) 2 + ( y ' - r' & - r 0 ‘ )*>(RI + 1 ) dz 
y xx ±2 


+ sin 


in*^ J A 


(x + r 6) dz + cos <p I A(y - r ©) dz 


I 


J- 


sin 2# I A(x + r S) ( y - r 0) dz 3 

y x 


<3. 1 ) 


where , 


I = A < 1-z) + i A ( l 2 -z 2 > + , 
l o 2 l 


A . , n+ i n+ 1 . 
n ( 1 -Z ) 


n+1 


I = 1 A (l 2 -z 2 > + i A (l S -z S ) 
z 2 o 3 t 


A . , n+2 n+2 . 

n ( 1 ~Z > 


n+2 


In the equation <3.1>, the -first two terms represent the bending 
and torsional energy terms respectively while the third term is the 
effect of blade rotation. 


The Reissner's functional is defined as C201 
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Reissner's functional is employed in order to eliminate the 
shortcomings of potential energy and complementary energy methods. 
While potential energy methods lead to good displacement field 
estimates and poor stress estimates, the complementary energy methods 
lead to good stress estimates but poor displacement estimates. 
Reissner's functional yields good stress and displacement fields 
simultaneously since it involves both the strain and complementary 
energy functions. 


I 
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The strain energy function for the blade is C193 
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The body forces are assumed to be negligible- Hence, 


JJLv* 


dv = 0 


(3.5) 


The nozzle excitation can be taken as bending forces F ( z , t ) and 
F^(z,t) in the x and y directions and twisting moment M(z,t). For a 
rotational frequency oi and n^number of nozzles the excitation is 
periodic with nozzle passing frequency v ( =n^o>) . Therefore, by 

Fourier expansion 
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The Reissner's dynamic functional is, therefore, 
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It is to be noted that all the above expressions the 
cross-sectional properties like A, I etc. are taken as polynomials of 
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which satisfy the boundary conditions 
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42 

v 

i 


h.h .1 
1 J 


r l j x y 

J 11 

J E< I 2 - I I ) 


dz 


x y x X y y 

1*1 11 11 


- f h.g; 


dz 


l h.h. 

J — H- 

o 


Q = 1 

l 

[ a f . 

dz 

Q = I 

l 

fa f . dz 

0 ^ ■ j 

l 

r b f . dz 

oi J 

1 Ox J 

n> , 1 J 

' mx J 

m+<S, 1 J 

p mx J 

Q = 1 

0 

1 

fa f . 

dz 

Q = I 

0 

1 

fa f . dz 

Q = j 

0 

1 

f b f . dz 

02 j 

o ° y J 

m, 2 J 

o 

m+tf , Z J 

o w,y J 

Q = 1 

l 

f a g . 

dz 

Q = 1 

i 

fa g . dz 

Q 

l 

f b g. dz 

03 j 

' OM y j 

m, 9 J 

J mM J 

m-KS,a 

• mx J 

o 



o 


O 


< 3 . 14 ) 


All other elements are zero. 
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3.2 DAMPING 

In the present study blade damping is modelled as a nonlinear 
-function o-f strain amplitude and speed of rotation for each mode of 
vibration. For ease of formulation, initially, damping is assumed to 
be viscous and proportional so that the energy dissipated is 

W = C T / p (3.15) 

d v 

The process of incorporating the nonlinear damping is considered 
later, in section 3*4. 

The equations of motion for forced damped vibration are 
CM3Cq> + CCKq) + CKKq) 

= <Q > + \ -CQ } cos mi>t + S fQ > sin mvt 

O Za m Zj mttf 

m m 

(3.16) 

where 

fq> T = < A. D . , . .F . > 

1 1 x 

The damping submatrices are 
C 

C = — M 

l i p n 

C 

C = — M etc. (3.17) 

12 p 12 

The rest of the elements are as given by equations (3.14). 

3.3 MODAL ANALYSIS 

Modal analysis is used to obtain the decoupled equations of 
forced motion. Using 

<q> - CUD Cr)> (3.18) 


the decoupled equations of motion are 
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CM3<r>> + cc3<r?> + ck3<t?> 

= <Q > +y <Q > cos mi>t > sin mi>t 

ON 4, mN 4, ( mttf > N 

m m 

<3. 19) 

where 


CM3 

= CU3 T 

CM3 

CU3 

CC3 

= CU3 T 

CC3 

EU3 

CK3 

- CU3 T 

CK3 

CU3 


and 


CQ > = CU3 T 

CQ > 

ON 

o 

f Q > = CU3 T 

f Q > 

mN 

m 

<Q > = CU3 T 

CQ 

< m+ > N 

m+ eS 


The individual uncoupled equations are of the form 


v** * 2 c k p kV p kV ■ 


Q 


ON 


2 


Q 


mN. 


cos mi>t + 


2 


mt «s > w 


sin mi->t 

(3.20) 


where 


K = K i / \ 

<k = c k / 2 vk 

for which the solution due to each term in the forcing function is 
written as 




Q / M, 
on. k 
k 


V. 


Q / M. 

mN. k 
k 


yAp* - <i»i>>*3* + 4C*p*(mv) J 


cos (mi^t — yj ) 


rn 
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Q / 

< m+ <5>N k 

k 


< nv t-<S> , 


sin (mvt - 


/ 


Cp‘ - (.mu) 2 ! 2 + 4cfp 2 <mu> 2 
*- k k 


v > 

m k 


where 


< 3 . 21 ) 


ys = tan 


2 W |W 

2 2 
p* - <mv> 


The response due to each harmonic can be written as 

Q / M 


H ■ 7) 
o 'o 
k k 


ON k 
k 


H -7) 
m k m k 


< tn*KS> , 


/ 


Q Z + Q 2 

mN <m-KJ>N, 

k k 


/ M 


y Cp 2 - <mv) 2 ] 2 + 4C*p*<mv) i 


cos ( mvt - y/ - 6 ) 


m k ro k 


< 3 . 22 ) 


where 



= tan 


-i 


k 


Q 

mN 

k 


The -following -formulae are used -for the computation of stresses 

M„ = M cos -rfr - M sin -tp- 
X x z y x 

M. . = M cos + M sin ifr 
Y y x x z 

The normal stress is 

M y X m x y 
°x = I + I 

yy xx 

The maximum shear stress occurs in the vicinity of one of the 
points where the largest inscribed circle touches the boundary of the 
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blade profile of the section C213. 


T 


r& 



C 


This stress is given by 


where 


K 



4 2 

1 + J F/AU 



in which 

du = elemental length along the camber line 
t = profile thickness as a function of u 
A = area of the section 


C = 


n 2 D 4 


1+ 0.15 < 


2 — 4 

n D 
16A 2 


D 

2r 


16 A 


if the radius of curvature of the profile at the point is 
positive 


C = 


2 — 4 

n D 


16 A' 




1+ C0. 118 In ( 1- 


2r 


> - 0.238 


D 


2r 


tan 


2<f> 

TT 

m 


if the radius of curvature of the profile at the point is 
negative 

D = diameter of the inscribed circle 

r = radius of curvature of the profile at the point 
<p = angle through which a tangent to the boundary rotates in 
turning or travelling around the reentrant position 

( radians) 

Stresses t and t are obtained from t using coordinate 

xx yx ra 
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transf or mat ion relations . 
The principal stresses are 


a = 0 

i 



(3.23) 


The basic mean stress (zeroth harmonic) and alternating stresses 
are determined from the principal stresses by using the formula 


a 

m 


/ 2 2 

C a + a - a a 
ait 2 3 2 3 


(3.24) 


3.4 DAMPING MODEL 

Damping is a critical parameter in the stress analysis of 
turbomachinery blading. It determines the stress levels expected in 
resonant vibrations. The main damping mechanisms influencing blade 
behaviour are interfacial damping at the root, material damping and 
gas / steam damping. 

Rao, Gupta and Vyas C123 conducted tests in vacuum spin rig and 
established the dependency of overall damping on the rotational speed 
and strain amplitude. In these tests, excitation was provided by 
electromagnets simulating nozzles in a turbomachine. From the 
filtered decay signals of each mode, equivalent modal damping ratios 
were obtained as function of rotational speed and strain amplitude. 
Typical equivalent modal damping values for the first mode obtained 
experimental ly by Vyas [193 are shown in Fig. 3. 8. 
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The earlier study by Vyas C191, however averaged the modal 
damping value over the range of vibration amplitude ( for the speed 
under consideration ) . He assumed that this average can be taken as a 
constant damping value which is valid for all amplitudes and employed 
it for the purpose of stress analysis. In the present study, an 
improvement has been made over the above approach by making a 
provision for the amplitude dependence of damping. The damping ratios 
for each mode can be expressed as polynomial functions of strain 
amplitude, e , at various rotor speeds. Thus, 

C = * o + € + a z + + a = e 5 (3.25) 

The constants a^ in the above equation can be obtained by a curve 
fitting routine. The strain amplitude, e, is however, itself 
determined by the damping, £ . Hence, an iterative procedure is 
adopted . 

Starting with an assumed value of strain and corresponding modal 
damping ratio, point 1, Fig. 3.2, the principal component of stress 
and resultant strain is estimated, point 1'. The average of the 
assumed and estimated values of strain and the corresponding modal 
damping ratio are taken as the starting values for the next 
iteration. The process is repeated till two successive values of 
stress are obtained within a specified level of accuracy. 

3.5 COMPUTER PROGRAM 

The computer program developed for studying the blade vibration 
is modular in form. Individual modules can be attached or removed 
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depending on the requirement. Each module can be run independently 
and in tandem with the main program. Starbase graphics library 
support has been taken to present the result in a graphical -form tor 
easier in terpretation . Fig. 3.3 gives the structure of the software 
developed. It is briefly outlined below. 

The package requires the following inputs to be fed by the 


user : 


1. Blade sectional profile data 

2. Blade material properties 

3. Blade damping data 

4. Nozzle excitation data 


■> 


-> 


•> 


•> 


sect.dat 


mat . dat 


damp_pl .dat 


excit .dat 


f lex.f 


reads data from fsect.dat] and computes the 
coordinates of the centre of flexure of the aerofoil 
cross-sec tion . 


sec t . f 


damp, f 


reads data from 


^ect .dat and ff lex .outj and determines 


the area, second area moments of the cross-sec tion 
about x and y axis, product moment of cross-section 
about x-y axis, polar moment of inertia about an 
axis passing through the centre of flexure and the 
torsional stiffness of the cross-section. 


reads standard values of damping ratio from 
|damp_pl .dat) and computes the damping curves for the 
particular speed at which the rotor is being run. It 
first locates that particular speed between two 
speeds for which damping data is known. It then 




FIG. 3.3 SOFTWARE 
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|dr_forc . f 


| stress . f 


interpolates the damping values and obtains the 
damping curve. 


reads excitation data from fexcit.dat and computes 
the Fourier force components of the forces and also 
plots them along with the forces. 


reads input from four files already generated, 
namely 

a) Section properties data f rom fsec t . out 

b) Excitation Fourier components from jexcit . in | 

c) Damping data for the rotor speed from [dp. out] 

d) Material property data from [stress, datj 

The computation proceeds along the following lines : 

1. Curve fits for the section properties with 
length is obtained. 

2. The mass and stiffness matrices are then 
generated using standard quadrature formulae. 

3. The eigenvalue matrix is then formulated. 

4. The eigenvalue matrix is solved and 
eigenvalues and eigenvectors determined. 

5. The natural frequencies and mode shapes are 
obtained using eigenvalues, eigenvectors and shape 
functions . 

6. The generalised stiffness, mass and force 
matrices are obtained. 

7. The decoupled equations of motion are then 
solved to obtain the response using any one value 
for the damping ratio ( the maximum in this 
implementation ) . 

8. The displacements are then determined and 
stresses obtained. The principal stresses are used 
to compute the strains, from which a new damping 
value is obtained ( since strain amplitudes and 
damping ratio are related nonlinearly ). If the old 
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and the new values of the damping ratios differ by a 
defined small quantity, then that is the correct 
value of damping and the stresses are reported. If 
they do not coincide, the average of the new and old 
values of the damping ratio is used and this step is 
again repeated. 


The files with the extension dr_ are essentially graphics generating 
program files. Those with extension tit_ are title files and contain 
the text to supplement the graphics. 


dr_prof . f 


reads profile data arranged in a particular format 
from [sec t . t r] and generates the blade profiles for 
various sections along the length. 


dr_sdamp. f 


reads standard damping data fed into jdamp _pl .dat [ and 
draws the standard damping data for each of the 
first four modes as a function of the rotor speeds. 


dr_damp . f 


dr modesh.f 


reads data about the damping curves for the 
particular speed at which the rotor is being run 
from [damp. out and draws the damping curves. 


reads data 
shapes. 


from Imodesh.outl and draws the mode 


dr strvsl.f 


dr_curv . f 

reads data from 

coef i .out 

and 

prop. out ] 


various section properties, mentioned earlier, as 
function of blade length. 


reads data f rom Istr .outl and draws the stresses as 


function of blade length for the various harmonics. 


tit _prof 


contains text for the drawing file jdr_prof . f 
title files operate similarly. 


Other 
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All the above programs are run by the file "start" which 
contains the commands for compiling and executing all the programs 
sequentially- Individual programs can be compiled individually and 
executed. New modules can also be added and/or deleted which is an 
important feature of the package. 

3.6 RESULTS AND DISCUSSION 

The package developed is run for an illustrative case for which 
data is taken from C193. Fig. 3. 4 shows the six cross-sectional 
profiles of a typical turbine blade at different distances from the 
root. The various geometric properties obtained from the profiles are 
plotted against blade length in Figs 3.5a and 3.5b. Typical nozzle 
excitation forces are periodic in nature, and as simulated 
experimental ly by means of electromagnets C193, are shown in Fig. 
3.6. The same figure also shows the Fourier component amplitudes of 
the forcing functions. It may be noted here that the magnitudes of 
the excitation forces are very small since they have been obtained 
from laboratory simulation by electromagnets. In practice, these 
forces are considerably large to drive the machine for power 
generation. The natural frequencies are listed in Table 3.1. The mode 
shapes are plotted in Figs. 3.7a through 3.7e. The first and second 
modes are predominantly flapwise bending modes, with the twist being 
extremely small. The third mode shows an increasing contribution due 
to twist. The fourth mode is predominantly torsional . The fifth mode 
is a coupled mode with deflection in the bending mode dominating. 
Figs. 3.8 through 3.11 show the damping curves for the first four 
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modes at different rotor speeds. Fig. 3.12 is the Campbell diagram 
for the blade with six nozzle passing harmonics interacting with the 
natural modes to give resonant rotor speeds upto 4000 rpm. The 
resonant speeds are listed in Table 3.2. Figs. 3.13 through 3.15 show 
the damping curves for the rotor speeds at which the fourth, fifth 
and sixth nozzle passing frequencies interact with the first mode. 
The blade stress fields have been computed for various rotor speeds 
upto 1200 rpm including resonances. The stresses are determined at 
six equidistant stations along the blade length. Fig. 3.16 shows the 
variation of the stress harmonics at the blade root as a function of 
the rotor speed. Figs. 3.17a and 3.17b depict the resonant stresses 
as a function of blade length. 

It can be observed that the stress fields on the blade are 
functions of the order of the nozzle passing harmonic, the amplitude 
of the forcing function components, the interacting blade natural 
mode and the blade geometry. For the resonant speed under 
consideration, ie. 930 rpm, the stresses are maximum at the root and 
decrease towards the tip. The maximum stress due to the interaction 
of the sixth harmonic with the first mode is greater than the 
corresponding value due to the interaction of the fifth harmonic with 
the first mode. This is because the Fourier amplitude of the fifth 
order harmonic is less than the sixth. A more rigorous parametric 
study can be undertaken, with the help of the computer package 
developed, for a detailed investigation of the operational and 
geometric parametric influences on the blade stresses. 
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3.7 CLOSURE 

In this chapter, the dynamic stress analysis o-f a tapered, 
twisted, asymmetrical aerofoil cross-section blade has been 
undertaken and the computer program developed illustrated for a 
typical blade under forced excitation and nonlinear damping. 





Fig. 3.5a Variation o-f geometrical properties with length 
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TABLE 3.1 


NATURAL FREQUENCIES <Hz) 

1. 743.2 

2. 2146.5 

3 . 4548 . 9 

4 . 7537 . 3 

5. 9743.2 

6. 13683.0 

7. 14475.0 

8 . 20884 . 0 

9 . 36682 . 0 
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Fig. 3.16 Stress harmonics 
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Fig. 3.17b Stress harmonics vs blade length at 930 rpm 


CHAPTER 4 


LIFE ESTIMATION 


Blade -fatigue is strongly influenced by the mean and alternating 
stresses experienced during operation, by the elastic and fatigue 
properties of the blade material and chemistry of the operating 
environment. The critical aspect of the problem constitutes 
occurrence of large dynamic stresses due to resonant or near resonant 
vibration at some engine harmonic. Large dynamic stresses, thus, get 
inflicted on the blade. The damaging fatigue influences cumulate to 
"initiate" a crack in a zone of high stress, at some metal lurgical or 
structural discontinuity, and if critical conditions of blade 
operations sustain, . the crack "propagates" to lead to failure. 
Reference can be made to article by Dewey and Rieger C223 for a 
survey of the nature and extent of some recent specific industrial 
turbine blade failures. 

It becomes imperative to correlate blade vibrations to blade 
fatigue in order to assess the residual life of the existing blading 
and for the development of newer blade designs. 

Prediction of blade residual life requires s 

(a) Stress loading history of the blade. 

(b) Material strength and life properties. 

(c) Life estimation algorithm. 
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These aspects have been discussed by Rao and Vyas C13 and blade 
life has been estimated using stress based approach ie. von Mises 
theory with the S-N mean stress diagram. Recently, Irretier 
presented a life estimation model based on different stress 
approaches to cumulative damage and discussed related numerical 
studies. 

These approaches, however, suffer from a drawback that they do 
not make allowance for the possibility of the development of plastic 
strain zones, especially in cases of low cycle fatigue. In this 
context, Rieger [173 and has pointed to the need for development of 
a fatigue model based on more fundamental and modern "fracture 
mechanics” concepts. 

In the present study, an attempt has been made to develop a life 
prediction algorithm based on fracture mechanics and strain-life 
approach. In order to duly analyse "crack initiation" and "crack 
propagation", a combination approach, originally proposed by Dowling 
[233 for notched members and modified by Socie [243, has been 
employed . 

The method consists of the following t 

Crack Initiation 

(1) Starting from a nominal stress ( from stress analysis ), 
determination of elasto-plastic strain using Nueber's rule 
£253 . 

(2) Coupling with material hysteresis curve and Massing 's 
hypothesis £263 and solving for plastic stress through 
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iterative technique. 

<3> Accounting for mean stress by Morrow's hypothesis [273 and 
determination of crack initiation life. 

Crack propagation 

(1) Determination of propagation life using Paris equation [283 . 

<2) Employing the method of Socie et al [243 in neglecting an 
estimate of the characteristic crack length. 

The method is outlined in more detail in Fig. 4.1. 

However, for the purpose of numerical demonstration of the 
method, a number of approximations have been made regarding the 
fatigue data. This has been essential due to lack of better 
criterion. The approximations have been stated explicitly whenever 
made. It is also to be realised that fatigue failure being a 
statistical process, experimental validation of life algorithms in 
real life situations proves to be prohibitively expensive. With the 
present state of art, real life turbine blading failure problems are 
being dealt only through approximate thumb rule predictions. 

In this context, the present model is to be viewed, not as an 
exact closed form solution of the problem but as an improvement on 
the existing conventional stress based approach. 

The various aspects of the methodology are described in the 


following sections. 
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4.1 STRAIN-LIFE APPROACH FOR CRACK INITIATION 

This approach is based on the observation that in many 
components, the response of the material in critical locations 
( notches and flaws ) is strain dependant. This is specially true for 
low cycle fatigue regime. In the strain life approach the plastic 
strain is quantified and therefore removes the drawback of the S-N 
approach which considers stresses and strains to be elastic. At long 
lives, where plastic strain is negligible, the stress life and 
strain- life approaches are essentially the same. 

Although most engineering components are designed such that 
nominal loads remain elastic, stress concentrations often cause 
plastic strains to develop in the vicinity of notches and cracks. Due 
to the constraint imposed by the elastically stressed material 
surrounding the plastic zone, the deformation at the notch is 
considered strain controlled. The strain-life method assumes that the 
smooth specimens tested under strain controlled conditions can 
simulate fatigue damage at the notch root of an engineering 
component. Crack growth is, however, not explicitly accounted for in 
this method and the need arises to go in for the fracture mechanics 
approach for crack propagation life estimates. 


4.1.1 Stress - strain relationships 

Total strain = Elastic strain + Plastic strain 

t + € (4.1) 

t e p 

For most metals, a log-log plot of the true stress versus true 
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plastic strain is modelled as a straight line. Therefore, 



and the expression for total strain is 



(4.2) 


(4.3) 


4.1.2 Cyclic stress - strain behaviour 

The response of a material subjected to cyclic inelastic loading 
is in the form of a hysteresis loop as shown in Fig. 4.2. The 
relation between the total strain range, elastic strain range and the 
plastic strain range is 

Ae = Ae + Ae 
e p 

or, 

Ae La + A<e 

_ P 

2 2E 2 

A comparison between the monotonic and cyclic stress— strain 
curve provides a quantitative assessment of the cyclically induced 
changes in the mechanical behaviour. The potential danger of using 
monotonic properties to predict cyclic strains is pointed out in 
[29 J . The monotonic properties may predict strains that are fully 
elastic, when in fact the material may experience large amounts of 
cyclic plastic strains. Therefore, whenever available, the cyclic 
stress— strain data should be used for cyclic loading. 


(4.4) 


(4.5) 
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4.1.3 Stress - plastic strain law 

Analogous to the monotonic stress-strain curve, a log-log plot 
of the completely reversed stabilised cyclic true stress versus true 
plastic strain can be approximated by a straight line and a similar 
power law function is used, 


O' = K ' ( e ) 
P 




(4.6) 


(4.7) 


The equation of the hysteresis loop can be derived from the 
equation of the cyclic stress-strain curve using Massing ' s 
hypothesis C263. This allows the stabilised hysteresis loop to be 
determined for a material that exhibits symmetric behaviour in 
tension and compression. This hypothesis states that the stabilised 
hysteresis loop may be obtained by doubling the cyclic stress— strain 
value from the stabilised cyclic stress-strain curve. The equation of 


the hysteresis loop can, now be written as 


[ 2K ' ) 


1/n ' 


(4.8) 


4.1.4 Geometry effects 

The fatigue resistance of most metals is affected by the 
geometric conditions. Stress concentration is evidenced due to 
surface discontinuities such as notches, holes and scratches. The 
stress concentration factor, K^_, is dependant on the geometry and the 
mode of loading and its values are compiled in various texts, the 
most popular and well used being that of Peterson's C303. 
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Tests of notched specimens usually result in fatigue lives which 
are greater than those using and for this reason a fatigue notch 
factor, K^, has been defined. In accordance with the approximations 
suggested by Socie C243, which will be discussed in a later section, 
the theoretical stress concentration factor, K^, will be used for 
formulation and computation. 


4.1.5 Nueber's rule 


Nueber C253 analysed a specific notch geometry and obtained his 
rule. Although his method was proved only for one notch geometry, it 
is assumed that his relationship holds true for most notch 

geometriesC29D . The rule states that the theoretical stress 
concentration is the geometric mean of the strain and stress 
concentration, ie. 

K. = ( K K > 1/2 <4. 9) 

t a e 

where , 

K ; stress concen tration factor = % 
a S 

K : strain concentration factor = 

<s e 


or, 

K? * K K <4.10) 

TL O € 

When nominal strains are elastic, the nominal stress, S, and 
nominal strain, e, are related by Hooke's law. The Nueber's rule then 


becomes 
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( K fc S ) 2 

■ a m < 4 . 11 ) 

I -f nominal strains are inelastic different version of Nueber's 
rule must be used. For the purpose of this analysis it is assumed that 
nominal strains are elastic. 



Basquin C31] observed that stress life < S-N ) data could be 
plotted linearly on a log-log scale. Hence, 

' < 2N f ) b (4.14) 

where , 

■ : true stress amplitude 

2N^. : reversals to failure 

o', ' t fatigue strength coefficient 

T 

b : fatigue strength exponent 

a ' and b are fatigue properties of the material , whose definition 
and determination is the subject of discussion of the immediate 


section . 
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Also, Coffin C321 and Manson C331, working independently, found 
that plastic strain-life data could be linearised on log-log scale. 
Therefore, 

Ae 

— = o' i ' ( 2N f ) <4.15) 

where. 


A«e 

— —■ : plastic strain amplitude 

' : fatigue ductility coefficient 

c : fatigue ductility exponent 

The total strain amplitude can now be written as 

| < 2N f > b + <= f ' < 2N f ) C 

The above equation (4.16) is known as Strain-life equation. 


<4.16) 


4.1.7 Determination of fatigue parameters 

The strain based approach to fatigue requires five 
parameters, namely a ^ ', b, c and E. 

Fatigue strength coefficient, ', is the intercept of the true 

stress amplitude - fatigue life < reversals to failure ) plot at one 

reversal . This is determined by curve fit of actual experimental 

data. In the absence of actual tests a fairly good approximation is 

a ' = a r < true stress at fracture ). For steels with hardness below 
f f 

500 HB it may be approximated as 

' £ S + 345 < in MPa ) 

f u 

where S is the ultimate tensile strength, 
u 
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Fatigue strength exponent, b, is the exponent obtained on curve 
fitting experimental cyclic fatigue data. Its value varies from -0.05 
to -0.12 for most metals with an average value of -0.085. 

Fatigue ductility coefficient, ', is approximated by the true 
fracture ductility ie. 

€ f “ € f = ln 1 - RA 

where RA is the percentage reduction in area. 

Fatigue ductility exponent, c, is not as well defined as others 
and a rule of thumb has to be followed for its estimation. Its value 
has been found by Coffin to be about -0.5, by Manson to be about -0.6 
and by Morrow to be lying between -0.5 and -0.7. The values of these 
parameters have been compiled such as those in C343. 

4.1.8 Mean stress effects 

Cyclic fatigue properties of a material are obtained from 
completely reversed , constant amplitude strain controlled tests. Mean 
stress is usually present in most practical situations. Mean stress 
can either increase the fatigue life, when compressive, or decrease 
it when tensile. 

The strain life equation is therefore modified, to include mean 
stress effects, as suggested by Morrow 1273 

A ^-- = — ■ ° * P . ( 2N ) b + e ' < 2N ) C (4.17) 

2 E p f P 

in which cr is the mean stress. 
m 

The predictions made with this equation are consistent with the 
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observation that mean stress effects are significant at low values of 
plastic strain, where elastic strains dominate. They also reflect the 
trend that mean stresses have little effect at shorter lives, where 
plastic strains are large. However, it incorrectly predicts that the 
ratio of elastic to plastic strain is dependant on the mean stress. 
Despite this drawback, it gives a fairly good life estimate. 

4.2 FRACTURE MECHANICS APPROACH FOR CRACK PROPAGATION 

Linear elastic fracture mechanics ( LEFM ) is based on the 
application of the theory of elasticity to bodies containing cracks 
or defects. As materials plastically deform, when the yield stress is 
exceeded, a plastic zone is formed near the crack tip. The basis of 
LEFM remains valid if this region of plasticity remains small in 
relation to the overall dimensions of the crack and the crack body. 

4.2.1 Stress Intensity Factor and Fracture Toughness 

Stress intensity factor defines the magnitude of the local 
stresses around a crack tip. THis factor depends on the loading, 
crack size, crack shape and geometric boundaries. Its general form is 
given by 

K = F /nST <4.18) 

where , 

a : remote stress applied to the component 

a : instantaneous crack length 

F s correction factor that depends on the specimen and crack 
geometry 
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Stress Intensity Factor solutions have been obtained -for a wide 
variety of problems and published in handbooks eg. C353. Finite 
element and other numerical methods are widely used for determination 
of K. 

As the stress intensity factor reaches a critical value, K , 

G 

unstable fracture occurs. This critical value of stress intensity 

factor, known as fracture toughness of the material, is a mechanical 

property measuring resistance to cracking under load. The plain 

strain fracture toughness, K , is dependant on specimen geometry and 

xc* 

metallurgical factors. ASTM E399, Standard Method of Test for Plain 
Strain Fracture Toughness of Metallic Materials, sets forth accepted 
procedures for determining this value. 


4.2.2 Fatigue crack growth 

The Paris equation C283 is used to predict the number of cycles 
spent in the growth of the crack to some specified length. In this 
equation 


C < AK ) 


(4.19) 


in which C and m are material constants and AK the stress intensity 
factor range. These material constants, C and m, are determined by 
curve fit of experimental data and are available for some steels in 
literature eg.C341. ASTM E647 sets guidelines for determining these. 


Using Paris formulation, the crack growth life, in terms of 


cycles to failure, may be written as 
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N 

P 



da 

C ( 6K ) m 


(4.20) 


i 

in which and a^ are the initial and -final crack lengths 
respectively. The equation (4.20) can be solved numerically. Before 
it can be solved the -final crack size, a^ , must be estimated using 



F varies with crack length, a, and iterative procedure are required 
to solve -for the final crack size. 


4.3 COMBINATION METHOD 

Dowling C233 proposed a method, which combines local strain 
approach to predict initiation life and fracture mechanics approach 
to predict propagation life, to estimate the total fatigue life of a 
notched component. He proposed that within a distance from the notch, 
1^, the local notch stress field dominates the stress intensity 
solution. When the crack is smaller than this length ( 1 < 1^> crack 
initiation can be estimated using the strain-life approach. Use of 
this approach to estimate the initiation life avoids the inherent 
difficulty of LEFM to describe short crack behaviour at the notch 
root. The basic assumption of LEFM, that the plastic zone size is 
small compared to the crack length, is violated for short cracks. 
Once the crack is larger than this length ( 1 > 1^_) crack growth can 
be modelled with the fracture mechanics approach. Therefore total 


life is given by 
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1 1 i © • Initiation life ( ) + Propagation life ( N ) 

<1 < l t > (1 > l t > 

STRAIN Lire APPROACH FRACTURE MECHANICS APPROACH 

(4.22) 

where the value of a^ for estimation of is given by 

a . = D + 1 . 
i n t 

in which D n is the depth of the notch . 

Dowling C233 has given estimates for 1^. 

A simpler combined approach that has been shown to yield very 
similar results to Dowling's method has been proposed by Socie C243. 
The implementation in the present work is based on Socie' s approach. 
His method combines 

• An initiation life estimate using the strain-life approach 
that assumes that the notch is fully effective ie. K^= K^. 

• A propagation life estimate obtained, using fracture mechanics 
approach, that assumes initial crack size equal to the depth of the 
notch. 

This method has simpler implementation since 1^ is not to be 

determined . 

4.4 IMPLEMENTATION OF THE ALGORITHM 

The algorithm for fatigue life estimation is depicted in the form 
of a flow chart in Fig.4.1. For the purpose of modelling, a small 
area in the vicinity of the root is considered as a plate. The width 
of the plate is taken as the maximum span of the aerofoil 
cross-section at the base. The plate is assumed to have the same 
section modulus as the aerofoil section and from this its thickness 
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is computed. It is assumed that the blade is defect free except for 
the surface discontinuity of specific nature as analysed in sample 
case studies. Such discontinuities will invariably occur during the 
manufacturing of the blade and are modelled as notches. The depth of 
such a notch is to be determined using non-destructive evaluation 
techniques for inspection of turbine blades and detection of cracks / 
flaws. If no crack is detected then the least count of the 
measurement technique can be taken as the initial depth of the notch. 
For illustrative case studies, a value of depth of notch equal to 1.0 
mm has been taken . 


From Fig. 3.16 which depicts the stress harmonics, it is seen 
that for rotor speeds upto 1200 rpm the root stress is maximum at 930 
rpm. It would, therefore, be of interest to estimate the life for the 
blade rotating at 930 rpm, at which the blade experiences resonant 
stresses due to the interaction between the first mode and the fourth 
harmonic of the excitation. The maximum section stresses occur at the 


root. The total mean stress on the blade comprises of the zeroth 
harmonic of the stress and the centrifugal stretching of the 


blade. The centrifugal stresses due to rotation can be determined as 


< R + z) 
A 


211 RPM 
60 


f ( J / d2 ] 


(4.23) 


A computer program based on the algorithm described is developed 
for estimation of the fatigue life. The blade material chosen is 
AISI-4340 steel and its properties ( from Ref.C14J > used f D r 
computation are listed in Table 4.1. The results for the case stucf j.es 
are presented in the next section. 
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TABLE 4.1 

PROPERTIES OF BLADE MATERIAL A I SI -4340 STEEL 
CYCLIC PROPERTIES 


Modulus of Elasticity ( E ) 

Cyclic Yield Strength ( S ) 

yc 

Cyclic Strain Hardening < n ' ) 
Exponent 

Fatigue Strength Exponent < b ) 
Fatigue Ductility Coefficient < ' ) 

Fatigue Ductility Exponent ( c ) 
MONOTONIC PROPERTIES 

Ultimate Tensile Strength ( S^ ) 

Yield Strength ( S y ) 

True Fracture stress ( a ^ ) 

Percentage Reduction < RA ) 

True fracture ductility ( ) 

Strain hardening Exponent < n ) 
Strength coefficient ( K ) 

MATERIAL CONSTANTS 

Fracture Toughness ( ) 

C 


193 GPa 
758 MPa 
0.14 


-0.076 

0.73 

-0.62 


1241 MPa 
1172 MPa 
1655 MPa 
57 7. 

0.84 
0.066 
1579 MPa 

137.375 MPa 
6.6 x 10~ 9 


m 


2.25 
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4.5 CASE STUDIES 


CASE i 

An initial crack, modelled as a circular notch o-f radius 1.0 
mm, is considered in the blade root section. This is depicted in Fig. 
4.3a. For 930 rpm,the stresses obtained -from the stress analysis are 
as -follows 

Alternating stress ( ) *= 0.056 MPa 

Mean stress due to zeroth harmonic ( a ) = 0.00966 MPa 

mv 

Mean stress due to centrifugal effects < a ) = 1.044 MPa 

me 

Total mean stress = cr + a — 1.05386 MPa 

mv me 

The width of the plate is determined as 3.8 cm and the thickness as 
1.0 cm. The theoretical stress concentration factor, K^., is 
determined as 3.0. The expression for correction factor F used, is 

F = [ sec [- 2 ~F-] ] 

The results of the blade life estimate obtained are as follows 

52 

Initiation life estimate = 1.7 x 10 cycles 

Propagation life estimate = 2.2 x 10 11 cycles 

52 

Total life estimate = 1.7 x 10 cycles 

The total life of the blade is very large since the stresses on 
the blade are small. This is due to the fact that the nozzle 
excitation forces due to electromagnetic simulation are very small. 
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In practice, the nozzle forces driving the machine are many times 
larger and inflict stresses which are considerably greater than those 
obtained, using available data, in the illustrated case. It is 
therefore assumed, for furthur illustration, that the forces are such 
that a thousand times larger stresses are inflicted on the blade. The 
stresses are now 

Alternating stress ( ) =56.0 MPa 

Mean stress due to zeroth harmonic (o' ) = 9.86 MPa 

mv 

Mean stress due to centrifugal effects < a ) = 1.044 MPa 

me 

Total mean stress = a + a = 10.904 MPa 

mv me 

The life estimates obtained are 

Initiation life estimate = 0.54 x 10*^ cycles 

Propagation life estimate = 0.32 x 10^ cycles 

Total life estimate = 0.54 x 10*^ cycles 

CASE 2 

Two initial cracks, modelled as two semicircular notches of 
depth 1.0 mm each, are considered in the blade root section. This is 
depicted in Fig. 4.3b. For 930 r pm, the stresses obtained from the 
stress analysis are the same as for Case 1. The theoretical stress 
concentration factor, K fc , is determined as 3.0. The expression for 
correction factor used is 

2 3 

F - 1.12 . 0.203(4] - 1.197 (4-] * 1.930 (-f-j 

The results of life estimate obtained are as follows 

13 v 

Initiation life estimate = 0 • 54 x 10 cycles 
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Propagation life estimate * 0.2 x 10 ^ cycles 

13 

Total life estimate = 0.54 x 10 cycles 

Comparison between Case 1 and Case 2 reveals that the crack 
initiation lives in the two cases are the same. This is because the 
phenomenon of crack initiation is critically dependent on the 
theoretical stress concentration factor, which is the same in the two 
cases since the ratio of flaw radius to the blade width is the same 
in the two cases. The propagation lives, though different, can be 
seen to be of the same order. On the basis of these observations, it 
can be stated the size of the flaw is more critical parameter than 
its location on the blade surface. Such an observation would suit a 
flaw detection technique whereby in case no flaw is detected, the 
least count of the detection technique is taken as the flaw size with 
its location being arbitrary. 


CASE 3 

An initial crack, modelled as a semicircular notch of depth 1.0 

mm, is considered in the blade root section. This is depicted in 

Fig. 4.3c. the expression for the correction factor used is 

2 3 4 

F * 1.12 - 0.231 [-£-] + 10.55 [- 5 -] - 21.72 (-J-] + 30.39 [-£-] 

The theoretical stress concentration factor, K t , is determined to be 
2.79. The results of life estimate obtained are as follows 
Initiation life estimate = 1.4 x 10 cycles 
Propagation life estimate =0.17 x 10 5 cycles 
Total life estimate = 1.4 x 10 13 cycles 
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CASE 4 

An initial elliptical surface crack into the three dimensional 
plate is considered (Fig. 4.3d). The correction factor for the stress 


intensity 

factor expression has 

been estimated 

using 

empirical 

relations 

from C351 . 

The stress 

intensity factor 

for 

elliptical 

cracks varies along the 

periphery 

and its maximum 

value 

has been 

taken for 

computation . 

Two cases, 

depending on the 

a/b 

ratio are 


considered . 

1) Semicircular surface crack with a/b ratio 1.0 

The values of a and b are taken as 1.0 mm each for illustration. 
The correction factor, F, is determined to be 0.732. The life 
estimates are 

13 

Initiation life estimate = 0.76 x 10 cycles 
Propagation life estimate = 0.47 x 10^ cycles 
Total life estimate = 0.76 x 10*^ cycles 

2 > Elliptical surface crack with a/b ratio 2.0 

The values for a and b are taken as 1.0 mm and 0.5 mm 
respectively. The value of the correction factor obtained is 0.657. 
The life estimates are 

Initiation life estimate = 0.65 x 10 cycles 

Propagation life estimate = 0.81 x 10 cycles 

12 

Total life estimate = 0.65 x 10 cycles 




Fig. 4.3e Conf iguration 


•for Case 5 ( not to scale ) 
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CASE 5 

An initial elliptical discontinuity is assumed to be embedded 
inside the plate considered as a three dimensional body (Fig.4.3e>. 
The value for the correction factor, F, is found from empirical 
relations given in C351. Two cases, based on the a/b ratio are 
considered . 

1) Elliptical embedded crack with a/b ratio 1.0 

The values for a and b are taken as 1.0 mm each. The maximum 
value of the correction factor is estimated as 0.638. The life 
estimates are 

13 

Initiation life estimate = 2.0 x 10 cycles 

Propagation life estimate = 5.9 x 10^ cycles 

Total life estimate = 2.03 x 10* 3 cycles 

2) Semi-circular embedded crack with a/b ratio 2.0 

The values of a and b are taken to be 1.0 mm and 0.5 mm 
respectively. The correction factor, F, is found to be 0.58. The life 

estimates are 

Initiation life estimate = 0.76 x 10 13 cycles 

Propagation life estimate = 0.78 x 10 3 cycles 

Total life estimate = 0.76 x 10* 3 cycles 
As in Cases 1 and 2, a comparison between Case 4 and Case 5 reveals 
that the order of the life estimates is the same irrespective of the 


location of the flaw 
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4.6 EFFECT OF VARIATION OF INITIAL FLAW SIZE ON LIFE ESTIMATES 

As seen in the earlier section that for a given loading 
condition, it is the flaw size that plays a crucial role in 
estimating the total life. A level of uncertainty can always be 
expected with any flaw detection technique while taking a value as 
initial flaw size for a new undamaged blade. With this in view the 
effect of initial flaw size has been studied. 

Case 3 in the previous section is adopted for i 1 lustration . The 
width of the plate is 3.8 cm and its thickness 1.0 cm. The mean and 
alternating stresses taken for this illustrative case are 99.6 MPa 
and 560.0 MPa respectively. The values for life estimates with crack 
sizes varying from 0.25 mm to 1.5 mm are tabulated in Table 4.2 . 

It is seen that the initiation life increases with 
increasing flaw size. This is due to the fact that with an 
increase in flaw radius there is a lowering of the stress 
concentration factor. Once initiated, a larger flaw can be 
seen to propagate more rapidly to failure. 

4.7 COMPARISON OF LIFE ESTIMATES BY THE PRESENT METHOD AND THE 
S-N APPROACH 

A comparison is made between the present method and the S-N 
approach < using the Bagci's failure surface ) as used by Rao and 
Vyas C 1 3 . 

4.7.1 BASCI'S FAT I SUE FAILURE SURFACE LINE 

Bagci's fatigue surface line is of the fourth order and is given 
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by 


°a-f " » 1 




(4.24) 


The S-N plane reads 


log S ' = - A log N + B 

S ' = S = k S 
e e u 


for 10 3 < N < 10 6 


for N > 10 


with 


,og 


(0.9 S >‘ 


B = log 


m r KJ r 

R < s~ ] 


k =* 0.5 for a standard rotating beam fatigue specimen made of 
steel . 

R^., the endurance limit modifying factor, is given by 


R = K K K K K K . 
f a b c d e f 


where 


K : surface factor 
a 

K. : size factor 
b 

K : reliability factor 
c 


K , : temperature factor 
d 

K • stress concentration factor 
e 

K, : miscellaneous effects factor 

T 


4.7.2 RESULT 

The blade material properties listed in Table 4.1 are used for 
comparison. It is assumed that a notch of radius 1.0 mm exists in the 
blade and the configuration corresponds to Fig. 4. 3c. The 
theoretical stress concentration factor for this geometry is found to 
be 2.79 ( from Peterson's curves ). The alternating stress and the 
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mean stress values used for the comparison are 85 MPa and 16 MPa 
respectively. 

The value of computed is 0.1 13, based on suitable values of 
the modifying factors. The life estimate using the Dagci's fatigue 
surface line is found to be 0.96 x 10 5 cycles. 

By the strain-life - fracture mechanics combination approach the 
life estimate obtained is 

10 

Initiation life estimate = 5.6 x 10 cycles 

Propagation life estimate = 6632.8 cycles 

Total life estimate = 5.6 x 10 l ® cycles 

It can be seen that the S-N approach gives a life lower than 
that obtained from the present combination approach, since it does 
not specifically account for the development of the plastic zone and 
directly assumes an initiated crack which propagates to failure. 

4.8 CLOSURE 

In this chapter, a combination approach to fatigue life 
estimation, based on the concepts of strain life and fracture 
mechanics, is outlined and its use illustrated for the estimation of 
fatigue life of turbine blades subject to alternating stresses due to 


nozzle excitation. 



CHAPTER 5 


CONCLUSIONS 


In this study the dynamic stress analysis of a typical tapered, 
twisted, asymmetrical aerofoil cross-section blade mounted on a disc 
rotating at a constant angular velocity has been undertaken. Stress 
analysis has been carried out for a blade with nonlinear damping 
properties. Damping has been defined as a function of mode of 
vibration, rotor speed and amplitude of vibration, and, a numerical 
procedure developed for incorporation of such damping in analysis has 
been presented. Results, obtained for illustrative cases, are 
presented as functions of rotor speed, excitation harmonics and blade 

length. 

A life prediction algorithm, duly incorporating crack initiation 
and crack propagation phenomenon, has been developed as an 

improvement on existing conventional stress based approaches, 
algorithm developed is based on the Combination Approach, whereby 
crack initiation life is appropriately estimated using Strain-Life 
Methodology while crack propagation life is estimated employi g 


Fracture Mechanics Concepts. 

, _ the l i-fe estimation 

The stress analysis code along wit 

algorithm forms the overall package. However, a comprehensive pac g 

would require as input the complete stress loading biography of the 
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blade. Blade histories would involve constant speed service 
operations along with variable speed operations as during the process 
of start-up and shut-down of the turbomachine. This study has 
restricted itself to constant speed rotor operations. Stress 
Estimation under nonlinear damping during variable speed operations 
like start-up/shut-down can be taken up as the next stage of 
analysis. Provision also needs to be made for incorporation of 
thermal stresses incurred by the blading during operation of the 
turbomachine . 

Life Prediction models are generally constrained in absence of 
experimental validation. Such validation is considered to be 
prohibitively expensive. However, there has been a strongly felt need 
to correlate already occurred blade failures to their loading 
histories. Collection of such real life blade data from user 
industries is essential to check upon the efficiency of a proposed 
failure prediction model. 
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